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Dirac electrons and domain walls: a realization in junctions of ferromagnets and 
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We study a system of Dirac electrons with finite density of charge carriers coupled to an external 
electromagnetic field in two spatial dimensions, with a domain wall (DW) mass term. The interface 
between a thin-film ferromagnet and a three-dimensional topological insulator provides a condensed- 
matter realization of this model, when an out-of-plane domain wall magnetization is coupled to 
the TI surface states. We show how, for films with very weak intrinsic in-plane anisotropies, the 
torque generated by the edge electronic current flowing along the DW competes with an effective 
in-plane anisotropy energy, induced by quantum fluctuations of the chiral electrons bound to the 
wall, in a mission to drive the internal angle of the DW from a Bloch configuration towards a Neel 
configuration. Both the edge current and the induced anisotropy contribute to stabilize the internal 
angle, so that for weak intrinsic in-plane anisotropies DW motion is still possible without suffering 
from an extremely early Walker breakdown. 

PACS numbers: 


I. INTRODUCTION 

Dirac fermions in two spatial dimensions have been the 
object of intense study in recent times in the condensed- 
matter world, e^ecially since the experimental realiza¬ 
tion of graphene^J'21 and, more recently, the discovery of 
three-dimensional topological insulators which 

host Dirac fermions as topologically protected surface 
states. On the other hand, magnetic DWs and their 
manipulation via applied currents and electromagnetic 
fields hold a prominent position in the field of spintron- 
ics, especially so since t he e xperimental realization of the 
“race-track” technologjl^. 

The search of efficient ways of moving a DW at the 
highest possible velocities has become of capital impor¬ 
tance. Manipulation based on the application of exter¬ 
nal magnetic fieldj^HUl^ current J^^^^^ and more recently 
magnonjSMlll and electric fieldj^SEI] has been proposed 
and, with the exception of magnonic manipulation, ex¬ 
perimentally realized. There is an upper limit on the DW 
velocity due to the phenomenon known as Walker break¬ 
down (WB^. Above a threshold applied current or mag¬ 
netic field, the internal structure of the DW, as described 
by its internal angle, becomes unstable. The net velocity 
of the DW is limited by this effect. As a consequence, the 
search of mechanisms that can stabilize th e inter nal an¬ 
gle of the DW has become an important taslP^EH. In this 
regard, stabilizati on m ediated by Dzyaloshinskii-Moriya 
(DM) interactionjSSmi and Rashba fieldj^ has been ex¬ 
plored and experimentally realized. 

The appearance of three-dimensional TIs has fo¬ 
cused the interest on what we could call Dirac-fermion- 
mediated ferromagnetisnP^ and spintronics. The as¬ 
pect of TIs that makes them highly valuable for spin- 
tronic applications is that the spin orientation of the 
surface electrons is fixed relative to their propagation 
direction, so that the effects of spin-orbit coupling are 
maximal. For this reason, some proposal have sug¬ 


gested that TIs could possess more efficient spin-orbit- 
induced torques than other materials previously consid¬ 
ered. Indeed, the strength of the spin-transfer torque 
per unit charge-current density, exerted by the TI surface 
states on the magnetization of an adjacent ferromagnetic 
permalloy thin film, has recently been measured to be 
great er tha n for any other spin-transfer source measured 
so fail^^EH 

In this context it seems natural and promising to 
study TIs coupled to ferromagnet and specifically to 
DWa^SHSll If couple a TI layer to an out-of-plane mag¬ 
netized ferromagnetic thin film containing a domain wall, 
the DW acts as a mass for the surface electrons. It is a 
realization of a system of Dirac fermions with a DW mass 
term. The theory of this system coupled to electromag- 
netism was studied in a field-theoretical context in Ref.l^Sl. 
These authors showed how the 2n-dimensional anomaly 
of the chiral fermions living in the DW is canceled by 
the anomaly due to the induced 2n + l topological mass 
term. This cancellation in three spacetime dimensions 
has been explicitly computed for fermions coupled to an 
abelian gauge field in Ref.l^, where the technical diffi¬ 
culties that appear when trying to obtain the effective 
action in the presence of a DW mass become apparent. 

In this Article, we obtain analytical expressions for the 
effective action for an external electromagnetic field of 
a system of Dirac fermions in two spatial dimensions, 
at a finite density of charge carriers and with a DW 
mass term. We look at a physical realization in junctions 
of three-dimensional TIs and ferromagnets with out-of- 
plane easy-axis anisotropy that host an out-of-plane DW. 
We show how the surface electrons of the TI induce 
an effective in-plane anisotropy energy, which stabilizes 
the DW in a Bloch configuration even in thin ferromag¬ 
netic films where in-plane intrinsic and dipolar (shape) 
anisotropies are relatively weak. Owing to the stabiliza¬ 
tion of the internal angle, DW motion is possible without 
suffering from a very early Walker breakdown. 
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We also show how equilibrium and nonequlibrium edge 
currents can be generated along the DW by applying, re¬ 
spectively, a gate voltage (doping with electrons/holes) or 
an electric field. This current exerts a torque on the mag¬ 
netization that drives the internal angle from the Bloch 
configuration towards a Neel configuration. We analyt¬ 
ically compute the behavior of the DW internal angle 
as a function of the chemical potential and find it to 
qualitatively agree with a recent numerical calculatiorf^. 
Furthermore, this edge current contributes to the stabi¬ 
lization of the internal degree of freedom of the DW in a 
similar way as the DM interaction. It turns out that the 
corresponding term in the effective action, although it is 
first order in the magnetization, qualitatively resembles 
that of an interfacial DM interaction. Here the DM in¬ 
teraction is tuneable, with the interaction strength given 
by the amount of edge current flowing along the DW. 

Our approach has the advantage that, since the calcu¬ 
lation is done at the microscopic level, the origin of the 
physics underlying the phenomenology can be traced and 
well understood. 


II. EFFECTIVE ACTION 

In this Section, we present a the calculation of the 
effective action for the electromagnetic field coupled to 
a system of Dirac electrons with a DW mass. Readers 
interested primarily in the physics and phenomenology 
of ferromagnet-TI junctions may skip to Sec. |III| 

Let us start from the action of 3d Dirac fermions (in 
the following nd will be used for n-dimensional) coupled 
to an external electromagnetic field: 

S = J (I) 

with the DW mass 

m = mo(7tanh(xi/5) (2) 

where x^ = (xq, Xi,X 2 ), xq is time, cr = ±1 and 6 are the 

topological charge and the width of the DW respectively, 
we define ^ = ^' 7 °, p. is the chemical potential, is 
given by 

/I 0 0\ 

= 0 u 0 (3) 

[o 0 vj 

and the gamma matrices satisfy the anticommutation re¬ 
lations { 7 ^, 7 “^} = 277 ^^. We work with a metric with 
signature (-I-, —, —, —) and with k = c = 1. We consider 
the general case where the velocity v is not necessarily 
equal to the speed of light c. 

We obtain the effective action for up to second 
order in the fields. The fermionic spectrum consists of 
a chiral massless state bound to the DW plus massive 
extended states (see Appendix 0 . If the DW is wide 


enough, massive bound states also appear. The total 
number of bound states is given by the largest integer 
less than A -I- I, with the parameter A given by 


(see Appendix [A|). Let us consider the case where we have 
a DW that is so steep that the only bound state is the 
chiral state (A < 1), so that we can do an enlightening 
separation. We can consider the system as described by 
two theories, one 2d edge theory, describing the chiral 
electrons localized near the DW center, and one 3d bulk 
theory describing the massive extended electrons. Each 
of these theories, considered in isolation, is anomalous. 
The 2d chiral edge theory, on the one hand, is well known 
to be chiral anomalouaSH, while the anomaly in the bulk 
theory is a consequence of the generation [via a Chern- 
Simons (CS) term] of a topological mass of opposite signs 
on either side of the DW. However, the anomalies cancel 
via the Callan-Harvey mechanisnP^, so that the com¬ 
plete theory is anomaly free. 

The perturbative calculation of the effective action 
would in principle require computing the fermionic prop¬ 
agator from the exact fermionic spectrum, and then per¬ 
forming the integration of the fermionic degrees of free¬ 
dom in the path integral. This was done to second order 
in and for 77 = 0 in Ref.l^ as an explicit verification 
of the anomaly cancellation, but these authors consid¬ 
ered only the case mo —>■ 00 and focused exclusively on 
those terms that contain either two-dimensional or three- 
dimensional antisymmetric tensors, which are the terms 
relevant for the anomalies. The complete analytic com¬ 
putation for finite mg and remains a formidable task. 

To make the calculation tractable, we must introduce 
a number of approximations. First, we consider the adia¬ 
batic limit, assuming a constant mass in the calculations 
and restoring the xi dependence at the end. This approx¬ 
imation captures the bulk contribution (extended states) 
and is reliable as long as the energy associated with the 
typical length of the inhomogeneities in the mass is much 
smaller than the energy of the extended states (^ mo). 
This means that the approximation is asymptotically ex¬ 
act, but near the DW center it translates into the con¬ 
dition v/S mo which is never fulfilled if we consider 
the case A < 1 . As a consequence, nonadiabatic correc¬ 
tions will appear near the DW. Furthermore, even if the 
condition is fulfilled, this approach cannot describe the 
contribution of the bound states. Hence, as a second ap¬ 
proximation, we add the contribution of the chiral state, 
assuming A < I to avoid further computations for the 
contribution of the massive bound states. As a third and 
final approximation, we compute non-adiabatic correc¬ 
tions to the CS term of the bulk contribution. To obtain 
these corrections, we impose gauge invariance and the 
cancellation of the anomalies in the two theories (bulk- 
edge correspondence). 
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A. Edge theory 

Let us first compute the 2d edge effective action. From 
the action of eq. 0 and the fermionic spectrum obtained 
in Appendix the classical action for the chiral mode 
can be written as: 

Sr,l = J d^x pl{xi)ip^^^^{xo,X2)(ido + eAo+ 


tadpole terms do not vanish. This gives linear terms in 
the effective action and an associated equilibrium (exter¬ 
nal electromagnetic fields set to zero) chiral edge current 
density along the DW. The linear terms can be computed 
from the tadpole-like term (in imaginary time): 


Teg — ie 


d^xp\{xi) {iAo + avA2 )x 


-I- av{id2 + 6 ^ 2 ) - p^^|;^^^^{xo,X2) 
where (see eq. ( A10[ )): 

Pxixi) = BQ{X)cosh^~^^{xi/6)x 


( 5 ) 


X 2^"! ^,A+i,i;-sinh^(a;i/(5) 


( 6 ) 


with Bo{X) defined in eq. (A15) and 2 F 1 the hypergeo¬ 
metric function. When tr = — 1 we have Sr and when 
CT = 1 we have Sr. 

Before proceeding, we partially fix the gauge. Let us 
set 


A^ 9^^'>{xi)A^{xo,X2) 


( 7 ) 


where is a given function of xi. The remaining 

gauge freedom is given by 


An —>■ An + dnU! 


( 8 ) 


with w 7^ U!{xi). 


1. Equilibrium edge current 

As a consequence of having a finite chemical potential, 
when integrating out the chiral fermion in eq. (0 the 


dqodq2 1 

dTT^ qQ + aivq2 — ip 


(9) 


By multiplying and dividing by qQ—aivq 2 +ip the integral 
in momenta can be rewritten as an even plus an odd part 
in qo- The integration of the odd part vanishes, while 
the integration over both momenta of the even part gives 
(going back to real time): 

Teq =J d^x pI{xi) iAo +( 7 VA 2 ) (10) 

SO that we can define the equilibrium edge current density 
as: 


jeq = “^ 7 'A(a;i) (l,cru) 


( 11 ) 


Here and from now on the Greek letters refer to dimen¬ 
sions xq, X 2 ' a = 0, 2. 


2. Chiral anomaly 

The action of eq. 0 is analogous to that of the 2d 
chiral Schwin ger model, which is well known to be chi¬ 
ral anomalous^, so that the gauge symmetry at the 
quantum level is broken. Integrating out the chiral 
fermionic degrees of freedom in eq. 0 up to second or¬ 
der in the electromagnetic field we get (using dimensional 
regularization)Pll 


anomaly 


Attv 


dxidxip\{xi)p\{x^)x 


dXodX2 Aa{xo, Xi,X2) T] - 


dagb 


2d^ 


cb ‘ 




r 


( 12 ) 


where da = {do,vd 2 ), Aa = {Ao,vA 2 ) and 5^ = r]‘^’’dadb. 
Note that chemical potential does not play a role here, 
and there are two reasons for this to happen. First for 
massless (non chiral) 2d fermions the theory at finite 
charge density is indistinguishable to that at zero den¬ 


sity. Second, the chiral anomaly is well known to be 
insensitive to chemical potential and temperatur^^^llll 
However, as we obtained in eq. (10) finite p plays a role 
at first order in for chiral 2d fermions. The effect is 
equivalent to applying a chiral chemical potential to non- 
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chiral and massless fermions, which activates the chiral 
magnetic effect in 2d generating an equilibrium current 
density analogous to that of eq. (11). 

Finally, from eq. (12) we can obtain the nonequilib¬ 
rium edge current density: 

^ , pI{x[)E2{xo,x[,X2) 


Jne = j dx\- 


do -I- avd2 


(l,crw) 


(13) 

with E 2 = doA 2 — ^2^0 the electric field in the X 2 direc¬ 
tion. 


B. Bulk theory 


the hrst terms in a derivative expansion, which will be 
useful later on when we treat the physical system. This 
expansion is justified in the low energy regime « rn?, 
breaking down when m? ^ p^. If < < TOq it turns out 
that this breakdown occurs in the region near the DW 
where the adiabatic limit is no longer reliable. Hence the 
validity of the derivative expansion coincides with the 
validity of the adiabatic approximation. 

Now let us look at the gauge variation of the full theory 
(edge plus bulk). The gauge variation of the edge theory 
(eq. ([^) is: 

dr = —— / d^xLup\{xi)E2 (19) 

47 r J 


To obtain the bulk contribution we will assume a con¬ 
stant mass for the fermions, restoring the xi dependence 
at the end of the calculations in what basically is an adi¬ 
abatic approximation, as we mentioned before. This way 
we have 3d Dirac fermions with a mass term that breaks 
time reversal symmetry. This system is well known to 
give a topological response under an external electromag¬ 
netic field in the form of a CS terirP^. 

To proceed we can always split the effective action into 
vacuum [p = 0) and matter (/i ^ 0) contributions, so 
that the matter contribution is zero at zero density: 


^bulk — Tq T ^matt 


(14) 


with: 

To,matt = 2 J (1®) 

where = {Aq,vAi,vA 2 ) and where 11^’' is the po¬ 
larization function. Doing this separation all ultraviolet 
divergences appear in the vacuum part, while the matter 
part remains finite. 


1. Vacuum contribution 


To be consistent with the calculations done for the edge 
theory, we will use dimensional regularization to treat 
the ultraviolet divergences. The computation of the 1- 
loop polarization function for p — 0 is straight forward. 
Separating it into it’s even and odd parts: 


Yrfiv _ yr/ai/ 
IIq — 11 q ^ 


n^o 


(16) 


and doing the computation we obtain: 
127ru^Vm^ 




d^d'' 


) (17) 


+ (18) 

with m = moCTtanh(a:i/(5). Here 9^ = {dcj,vdi,vd 2 ) and 
q 2 _ Note that we presented the results as 


which should be canceled by the gauge variation of the 
CS term (eq. (dH)) so that the theory is gauge invariant. 
This is not the case if E 2 is a function of Xi. The rea¬ 
son is that while the nature of the edge theory is totally 
non-adiabatic, the bulk theory has been computed in the 
adiabatic approximation. To cure this issue, we write the 
CS term as: 


/dlV 

0, 


e^aF^(xi) 

47rv^ 


e^^P'^dp 


( 20 ) 


so that the non-adiabatic term F\{xi) is fixed by impos¬ 
ing the anomaly cancellation. This is done in Appendix 
where the explicit form of F\{xi) can be found. 


2. Matter contribution 


The computation of the matter contribution is more 
involved. We again do the separation: 


-1- 

(21) 

The even part is computed in Appendix where adding 
both vacuum and matter contributions we get: 

\2ii\m\ 


+ e{P 

(22) 

rrOi _ rriO _ af^2 ,,2^ e a 

— -^-^e — y(jTl /i j . . 

127vv\m\ 

(23) 


(24) 


Some approximations have been done to arrive to these 
expressions. First, we are in both the adiabatic and 
the low energy limits. Second, the result is obtained by 
adding the polarization functions computed in both the 
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static {po — )■ 0) and the long wavelength (p —>• 0) lim¬ 
its. As explained in Appendix [Cl this means that within 
our approximations non-local terms which are constant 
in the limit po 0 and zero in the limit p —>■ 0 are being 
approximated by the constant term in 11°*’. Third, we as¬ 
sume that in the static limit the spatial momentum |p| is 
smaller than the Fermi momentum pp. This is generally 
true at low energies, except when the Fermi energy |/r| 
is above but very close to the value of |mo|j so that pp 
is very small. In this situation some corrections which 
are highly non-local would contribute. And fourth, we 
neglect dynamical contributions which are non-local (in¬ 
verse powers in the spatial momentum), which we believe 
will not have an appreciable effect in the description of 
the physical system, as we acknowledge in Appendix |D 2| 
To complete the bulk part of the effective action let 
us turn our attention to the odd part of the polariza¬ 
tion function, this is to the matter correction to the CS 
term. The calculation is straight forward and we obtain 
(vacuum plus matter contributions): 

nr = 


+ (25) 


rest of the main text): 


SpM = d 


—(^(cos6» - 1) - Hpm) 

7 ^ 


(26) 


which is the sum of the Berry phase term (“kinetic energy 
of spin precession”) plus the Hamiltonian: 


1 

\J 

9m 

2 

9m 

. J 

JdpM = 2 


dx 

+ 

dy 

j — Kml -1- Kpml 


(27) 

(see for exampl^^UZol) Here d is the film thickness, 
7 = pBge/fi the gyromagnetic ratio {ge = 2), Ms the 
saturation magnetization, A the exchange constant (ex¬ 
change energy per unit length) and K and Kp the easy 
axis and hard axis anisotropy constants (anisotropy en¬ 
ergy per unit volume). The magnetization unit vector is 
m = (sin0cos((), sin0sin((),cosd). It relates to the total 
magnetization as M = hy m/a° = Mg m (where a is the 
lattice constant) and couples to the spin of the surface 
electrons of the TI insulator via an exchange interaction. 

The action for the TI surface electrons, including the 
coupling to the magnetization, takes the standard form: 


with 


Sti 


J dtdxdy — Hpi — 


(28) 


Note that all the non-adiabaticity is encoded in F\{xi), 
whereas no non-adiabatic corrections have been com¬ 
puted for p? > w?. The reason is that corrections in this 
last case can not be computed as we did for p? < as 
the anomaly cancellation can not be invoked. The gauge 
non-invariance of the CS term for p^ > can not be 
canceled by the chiral anomaly, which is i nsens itive to the 
chemical potential (and temperature T)l °2|63| Further¬ 
more, for finite p and/or T there appear infinitely many 
terms in the perturbative series (in A^) that break gauge 
invariance in the presence of a boundary (or a DW mass 
as in our case). These terms are parity breaking (the CS 
is the lowest order of this terms), and in the presence of a 
boundary are gauge invariant for zero p and T (the only 
exception to this is the CS term) and gauge non-invariant 
for finite p and/or T (se^^for a computation of the next 
order parity breaking term for finite T). Therefore gauge 
invariance can only be restored at the non-perturbative 
level when all the terms in the perturbative series are 
summed up, in a similar way as occurs with large gauge 
invariance in the theory with no boundarie^SMSfi] 


III. JUNCTION OF A FERROMAGNET AND A 
TOPOLOGICAL INSULATOR 

Let us now look at a condensed-matter realization of 
the theory of Sec [H] We will take a thin film of a fer- 
romagnet and place it on top of a three-dimensional TI. 
The action for the ferromagnet is (we restore h for the 


Hpi = vpz ■ {i{hV — ieA) x (t)± 




'xy 


± A^m^az - eAo, 


(29) 


where the surface of the TI is taken to be in the z = 0 
plane, and where we have defined the density to be zero 
at half-filling. Aj:y and Aj, are the in-plane and out- 
of-plane exchange couplings respectively (both are defi¬ 
nite positive); vp is the Fermi velocity of the electrons; 
m = {m^,my,mz) = {m.a;y,mz) and cr = {(Jx,cFy,cFz) = 
{(Ts;y,<Jz)] Ay = (Ao,A) and the signs -I- and — corre¬ 
spond to antiferromagnetic and ferromagnetic exchange 
couplings, respectively. Spi can be written as eq. Q 
(with the representation for the Gamma matrices given 
by (|A2|)) setting x^ = {t,x,y), v = vp, m = 


and doing the substitution: eAy ea^ = (eAo , eA^ ± 
AxyXriy j Vp , eAy T AxyTx ix I f') • 

The Hamiltonian (271 supports a non trivial minimum 


energy configuration in the form of the well known Bloch 
DW: 


00 = 2 arctan(e'^'^/‘^), (/o = ±7r/2 (30) 

where 6 = ^AjK and cr = ±1 are the width and the 
topological charge of the DW respectively. This is valid 
if K± 7^ 0, while if Kp = 0 the vacuum would be degen¬ 
erate and (j)o could take any value. This configuration 
gives the equilibrium magnetization: 

= (^sech{x/S) cos(j)o,sech(x/6)sm(j)o,ataTih(x/d)) 

(31) 
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We will work with fluctuations around the equilibrium 
configuration: 6{t,x) = 9o{x) +6{t,x) and (l){t,x) = (j)o + 
(j){t,x), giving the total magnetization: 

m = (sin 0 cost/), sin d sin (/),(Ttanh(a:/(5) + fhz) (32) 

with rhj = cos 6 — cos 9q. Note that we have imposed the 
magnetization to be homogeneous in the y direction. 

The gauge fixing condition given by eq. 0 has now 
to be fulfilled by the effective vector potential a^. We fix 
the gauge: 


eoo = —eExX — eEyy 

(33) 

eax = ±Axymy{t,x)/vF 

(34) 

eay = TAxymx(t,x)/vF 

(35) 

(t, x) = sin(do(a:) + ^(i)) cos (/i(t) 

(36) 

(t, x) = sin(0o(a^) + ^(i)) sin(/)(t) 

(37) 


where we have chosen an electrostatic configuration of 
the electromagnetic field {E^^y are constants). Since the 
X dependence of nixy has been fixed, the fluctuations can 
only be functions of t. To obtain an effective action for x 
dependent fluctuations we have to relay on the adiabatic 


approximation, restoring the x dependence at the end. 
For the adiabatic approximation to be valid the wave¬ 
length of the spin waves Isw in ths x direction has to 
be much bigger than the typical wavelength of the sur¬ 
face electrons 4; = fiyp/Az- Assuming the parameter 
A = AzS/{Hvf) introduced in eq. Q that defines the 
number of existent bound surface states (see Appendix 
0 to be A ~ 1, we have Isw >> S. 

This way, we have ’’almost” a completely analogous 
theory for the fermionic sector to that in the previous 
section, but now for the effective electromagnetic field 
a^. We say ’’almost” because in addition we have an 
extra field with which we have to deal. Before pro¬ 
ceeding, let us fix the values of the parameters. We set 
(the MI parameters are obtained froirPi^): d = 3nm, 
Ms = 3x lOM/m, A = 10"“J/m, K = 2x lO^J/m^, 
vp = 5 X lO^m/s, Az = Axy = 30meV. So that: 
6 « 7.07nm, A = 0.672. We will assume a very small 
perpendicular anisotropy K±, so that it can be neglected 
compared to the effective anisotropy induced by the TI 
surface electrons (see next two paragraphs and Appendix 
0. 

The computation of the effective action for the mag¬ 
netization is done in Appendixes 01 relaying on the 
calculations of Sec. im The total action reads: 

r = Sfm + T'ti (38) 


with: 


Tti 


j dtdxdy 


Axy S 


m 


xy 


- y J dx' ruxit, x)p\{x)pi{x')mx{t, x') 


(39) 


and where the spin density s is: 


s 




2hvF 


aFx{x) d(A^ tanh^(a;/i5) — /i^) 


aAz tan\i{x/5) 

IaI 


6»(Ai^-A^tanh^(x/(5)) E 


(40) 


Here AV[y) = V{y) — V{—Ll2) is the voltage between 
y = —L/2 and a given point y along the DW, with 
r(y) = —Eyy, and L and y = ±L/2 the length and the 
end points of the wall respectively. The effective hard 
axis anisotropy constant is: 

Reff ^ ^ (41) 

dohvF 

and the functions px{x) and F\{x) are given by eqs. 0 


and ( B4|B5 1 respectively. Remember that the parameter 
A is given by (see eq. 0): 


A = 


A.d 

hvF 


(42) 


It is important to note also that the spin density is related 
to the electromagnetic current density as j = ieui^’S x z. 

From eqs. ( 39|40 1 we see that there is a spin density 
that couples to the magnetization, which is related to: 1) 
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(term proportional to x in eq. ( |40[ )) the edge equilibrium 
and nonequilibrium currents flowing along the DW, given 
by eqs. (11) and ( |D5[ ) respectively and 2) (term propor¬ 
tional to E in eq. (|40p) the topological current due to the 
anomalous quantum Hall effect in the bulk, coming from 
the Chern-Simons term (251. Besides the spin density 


there is a non-local contribution induced by the chiral 
electrons bound to the DW, which acts as a hard axis ef¬ 
fective anisotropy energy in the direction perpendicular 
to the wall. 

Although a interfacial DM interaction term of the form 
m • V TTiz is not induced, as argued in Appendix the 
spin density related to the edge current generates a term 
that resembles it: 




(43) 


For more clarity, in the special case A = 1 we can write 
this as: 


y - eAV{y) 
2hvF 


mz;y ■ 


(44) 


while for different values of A small deviations from 
take place. The strength of this ’’pseudo” DM 
interaction can be tuned by doping with electrons/holes 
and by the aplication of a voltage between both end 
points of the DW. Note however that this term is first 
order in the magnetization, while the DM interaction is 
second order. 

In Fig. [2 we show the schematics of a possible experi¬ 
mental setup. A thin ferromagnetic film hosting a DW is 
deposited on top of a TI, which itself is deposited on top 
of a gate. Two electrodes are attached to the flanks of 
the TI such that a voltage between both edges is applied, 
and the gate is used to tune the chemical potential in the 
TI. A similar setup was realized irP^ where the TI was 
used as the channel of a field effect transistor. 


IV. PHENOMENOLOGICAL RESULTS 



FIG. 1: Schematics of a possible experimental setup. The 
ferromagnetic thin film (green) is deposited on top of the TI 
(blue) which itself is deposited on top of a gate. Pink and 
blue thick arrows represent the out-of-plane magnetization 
pointing up and down respectively. A voltage AV is measured 
between the electrodes (yellow) at the flanks of the TI. The 
edge current ledge flowing along the DW is represented by the 
black arrow connecting both electrodes. 


From the spin density of eq. (401 we see that an applied 


electric field in the x direction will not change the average 
chirality of the wall, since it contributes with terms anti¬ 
symmetric in X which will vanish when integrated. Con- 
trarily, a voltage in the y direction has a non vanishing 
effect on the average chirality through the (nonequilib¬ 
rium) edge current that is generated (term proportional 
to the voltage in eq. (40)). This adds up to the effect of 


doping the system with electrons/holes, which generates 
a further (equilibrium) contribution to the edge current 
(term proportional to the chemical potential in eq. (40)). 

Then at finite density and under an applied voltage 
in the y direction we obtain the potential energy for the 
DW (after integration in x and y): 


E{<P) = 

hVF 


— Cl cos (j) + 


AxyC2 

2 


cos^ 4 > 


(45) 


with: 


Now that we have the effective action we can study 
the phenomenology. We will compute the chirality of the 
wall in it’s equilibrium configuration at finite density and 
under applied external electric fields. We will also look at 
the current induced dynamics when a current is applied 
through the ferromagnet. 


A. Equilibrium configuration of the DW 


We want to look at minimum energy configurations of 
the total effective action given by eqs. ( 38|39| ) of the type: 
m = {sech{x/5) cos (j), sech{x/5) simp, crtanh(a:/<5)). In 
principle (/ will be a function of x, however to simplify 
things we will assume (j) to be constant. This will give an 
average equilibrium value of 4‘{x). 


Cl = ay, — 


eAV 


(46) 


C. = -2i 


6_i(I/2 + A,-2A) 
e_i(A,l-2A) 


(47) 


where AV = V{Ll2) — U(—L/2) is the voltage between 
both end points of the DW, B is the incomplete Beta 
function (see eq. ( A16| )) and A is given in eq. (42). Here 
Cl can take negative values while C 2 is real and always 
positive and fulfills lim;),_,,oo C'2 = 1. For |Ci| < AxyC 2 
the energy has a minimum at ^ = arccos(C'i/(Aa;j,C'2)). 
On the other hand, for jCil > AxyC 2 the minimum en¬ 
ergy configuration is ^ = 0 for Ci >0 and cj) = n for 
Cl < 0. So at Cl = 0 we have a Bloch DW (</ = ±7r/2), 
and as |Ci| increases (j) is shifted until |Ci| = AxyC 2 , at 























( 48 ) 



LC2^xyIedge 2 j, j. tp v 

- — COS 4> - Teicj) - FelX 

evp 

where Af = 2d6L/aA is the number of spins in the DW, 
the currents ledge and I^^ge are: 

^ e , elsy. 

ledge = ^ (ct/T ^ ) (49) 

pA 

ndge = C2^^^-l0xlQ-^A (50) 

and Tel and Fei are the spin-transfer torque and t he for ce 
generated by the current on the wall, respectivelj^^^IM 


FIG. 2: Equilibrium value of the internal angle as a func¬ 
tion of the chemical potential and the voltage, for a A param¬ 
eter value of A = 0.672. 


which point the DW estabilizes in a Neel configuration, 
with (/) = 0 for positive Ci and (/) = tt for negative Ci (see 
Fig. §. This way the chirality of the DW can be tuned by 
the chemical potential (applying a gate voltage) and the 
electric field (appliying a voltage between both end points 
of the DW). Based on this mechanism, there is a way to 
switch between the two degenerate vacua (f) = ±7r/2 using 
an out-of-plane magnetic field to break the degeneracy of 
the vacuum. It was described iiP, where they made use 
of an electric field to tune the chirality through the mech¬ 
anism explained above, while here we have shown that it 
can be done also by applying a gate voltage. 

To be specific let us choose the topological charge to 
be cr = -1-1 and switch off the electric field. For the 
parameters introduced in sec. |T^ A takes the value 
A = 0.672, so that we have: Ci = fi, C 2 ^ 0.706. For 
1^1 < 0.706A3 ;j^ we get 4> = arccos(/r/(0.706Aa;y)) while 
for 1^1 > 0.706Aa,y we get <() = 0 for positive /r (TI doped 
with electrons) and (j) = n for negative /i (TI doped with 
holes). This behavior qualitatively reproduces the results 
obtained numerically irPSl. The physical explanation lies 
in the competition of the torque generated by the chiral 
edge current on the magnetization and the effective hard 
axis anisotropy energy, which is also induced by the TI 
chiral surface electrons. 


B. Current induced magnetization dynamics 


Tel = ^Is ; Fel = -JIs (51) 

e eo 

Here Is = — Ii is the applied spin current through the 

ferromagnet in the positive x direction, and the j3 term 
is a constant that depends on the microscopic properties 
of the ferromagnet. It can be quite large for perpendic¬ 
ular anisotropy thin films, where the force from electron 
reflection can be dominanlP^, so following the reference 
we will fix it to be /3 = 0.3. 

Now we can write the equations of motion (we include 
Gilbert damping): 

X LC 

eM{—- a(j)) = {ledge sin (j) - sii\2(t>)-Is (52) 


eN{4'+ a^) =-P Is (53) 

where a is the Gilbert damping parameter, which we will 
set to be a = 0.01. From the previous expressions we can 
obtain a differential equation for (j): 


where Is 


Js = 


3 s 

= dLjs, It = 
ct' = 


eSCi e,f.f(l3 
2h 


-I- sin 2(j) 

^edqe • j 

„ sm<() 

(54) 


'^edge 


dLjt and: 


1 + 


(55) 



-1)"' 

« 3.39 X 10^A/m^ 

(56) 


Let us apply a current through the ferromagnet in the 
direction perpendicular to the DW. In terms of the col¬ 
lective coordinates X{t) and (j){t), where X defines the 
position of the DW, we can write the total effective ac¬ 
tion as: 

V [ J4.(A I djC2XxyIedge , 

F = / dt\Mh—(j)-\ - - - —coscj)— 

J \ 0 evF 


For vanishing edge current ledge (m = 0 and AH = 0) 
the DW moves with a time-averaged terminal velocity 
{X) = —I361s/{eaj\f) as long as the current flowing 
through the ferromagnet is smaller than the critical cur¬ 
rent /*. When Is reaches the value /* however, Walker 
breakdown (WB) occurs and (j) starts to change, so that 
the averaged terminal velocity decreases as Is increases 
(see for exampl^^. Thus to obtain high velocities it is 
important to stay in the non WB regime, which means 
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FIG. 3: Density plot of {<j}) in the ledge/Ildg^ ~ js/j/ plane. 
Values {(j)) < 0 (WB) occur in the dark region, while in the 
light region {</) = 0 (no WB). The frontier between the two 
regions gives the maximum values of js before WB occurs. 


that one should look for the biggest possible I*. We see 
that it’s value is proportional to so that the big¬ 
ger the exchange coupling the higher the velocities that 
can be achieved. For the actual values of the parame¬ 
ters a velocity {X) « 1.86m/s is achieved for an applied 
current density js = jt ■ However if one could achieve an 
increase of A^y of one order of magnitude {Ag;y ~ 0.3eH) 
the average maximum velocity would increase by two or¬ 
ders of magnitude: {X) « 186m/s. This way we see 
that by using a TI, ferromagnets with very week hard 
axis anisotropy are still suitable for hosting DW motion 
without suffering from early WB. 

There is a nice way to stabilize the chirality of the DW 
so that WB is delayed. It was first suggested irP^ where 
an electric field was needed. The chirality stabilization 
is mediated by the edge current through the third term 
of the r.h.s. of eq. (54|. We plot in Fig. ^ the val¬ 
ues of js at which WB occurs for different values of the 
edge current ledge- Notice that DW stabilization can be 
achieved even in the absence of an electric field by doping 
the TI with electrons or holes through the application of 
a gate voltage. Indeed both mechanisms can be added so 
that higher edge currents can be achieved, and therefore 
higher terminal velocities can be reached. 


V. CONCLUSIONS 

We have analytically computed an approximate effec¬ 
tive action up to second order in the electromagnetic 
field, for a system of Dirac electrons in two spatial di¬ 
mensions at finite density and with a DW mass term. 
We have presented a condensed-matter realization of this 
system, consisting on a ferromagnet hosting a DW (with 
out-of-plane magnetization domains) coupled via an ex¬ 
change interaction to the TI surface electrons. There are 
three relevant contributions to the effective action com¬ 
ing from quantum fluctuations of the fermionic surface 


states: 

• The first one is linear in the magnetization, which 
couples to a spin density in the perpendicular di¬ 
rection to the DW. This spin density is related to 
the edge current of chiral electrons flowing along 
the DW, which itself can be seen as made of two 
pieces: an equilibrium current proportional to the 
chemical potential and a nonequilibrium current 
proportional to the applied voltage between both 
end points of the DW. 

• The second contribution is again linear in the mag¬ 
netization, but now the spin density to which it is 
coupled is proportional to and has the direction of 
the applied electric field, and is of opposite signs to 
either side of the DW. It is related to the topolog¬ 
ical current generated by the anomalous quantum 
Hall effect in the bulk. 

• The last contribution is a non-local term quadratic 
in the magnetization, induced by the chiral elec¬ 
trons, and which acts as an effective hard axis 
anisotropy energy in the direction perpendicular to 
the DW. 

The competition of the torque exerted by the edge 
current (first contribution) and the effective hard axis 
anisotropy energy (third contribution) explains the be¬ 
havior of the chirality of the DW as the chemical poten¬ 
tial and/or the voltage between both end points of the 
wall are modified. The stabilization of the internal angle 
through the induced effective hard axis anisotropy allows 
for the motion of the DW with velocity proportional to 
the applied current through the ferromagnet, especially 
in the case of ferromagnetic thin films with very weak 
in-plane anisotropy, which other ways would suffer of a 
very early WB. The critical current at which WB occurs 
is proportional to the effective anisotropy energy, which 
itself is quadratic in the exchange coupling, so that in¬ 
creasing the value of the exchange coupling would result 
in a significant increase of the maximum possible DW 
velocity. 

Finally, the edge current flowing along the DW has an 
interesting effect on the wall dynamics. It further sta¬ 
bilizes the internal angle of the DW, which translates in 
a delay of the appearance of WB. This means the max¬ 
imum DW velocity can be increased by applying a gate 
voltage (doping with electrons or holes) and/or apply¬ 
ing a voltage between both end points of the DW. The 
two effects can be combined so that the edge current is 
further increased and the WB further delayed. 
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Appendix A: Fermionic spectrum 

To obtain the fermionic spectrum, we will treat the 
gauge field as a perturbation and think of the remaining 
theory (eq. Q) as a free theory. For the free theory, 
the equations of motion can be obtained from eq. Q for 
vanishing electromagnetic field: 




(Al) 


Let us take the following representation for the gamma 
matrices: 


.,,0 ^3 

j = a , 


y.2=_,^1.2 


(A2) 


and write the bispinor in the basis: 

(A3) 


with: 


1 A 


V2 VI 


ul = 


A2 1-1 


(A4) 


so that -I- With this representation and 

in this basis, rewriting the first order coupled equations 
(Al) as second order decoupled ones we get: 


(a" + sech(xi/S)^ = (^ +1^2 - ^) 

(A5a) 


sech(xi/S))<l>L = (^+pI- 


■ _i_ 

g2 

(A5b) 

with A = tuqS/v. This is nothing but the Schrodinger 
equation in a modified Poschl-Teller potential (see for 
exam pl^. Let us writ e Pn = 

the eigenvalues of eqs. ( A5a | A5bt are —k^. The general 
solutions of these equations ar4^ 


$fl(cr,a:i) = cosh^+ 2 {xi/5){B{\)2Fi 


1 — a 


, b + 


1 — a 1 


; — sinA(a;i/<5) 


+ iC'(A) sinh(a;i/(5) 2 F 1 


1 — Cr 1, 1 — (T 13 .,2/ 

a+ + ^,b+ -;-smh {xi/5) 


r 


(A6) 


and ^Licr^Xi) = ^R^—a^Xi). Here B{X) and C'(A) are 
arbitrary constants, 2 F 1 is the hypergeometric function, 
and the values of a and b are: 


for the even part of eq. (A6) and: 


a=--{X + i6k), b=-{X — iSk) (A7) 


Note that the first term of the r.h.s. of eq. (A6) is even 
in xi while the second term is odd. For > 0 we have 
the continuum of extended states, which general solution 
is precisely given by eq. (A6). On the other hand, for 
k^ < 0 the solutions are bound to the wall, vanishing at 
±ocP, and the energies are quantized. 

Let us have a close look at the bound states. If we 
set k ^ ik we obtain the following conditions for the 
solutions to be normalizable (for a = —1): 


—>■ Sk = X — 2n 
—> 5k = X — 1 — 2n 


(AS) 


q>L 


5k = X — 1 — 2n 
5k = X — 2 — 2n 


(A9) 


for the odd part, where n = 0,1, 2,... (for tr = 1 one just 
has to interchange the chiralities: ^r^l ^l,r)- From 
these normalizability conditions and the general solution 
of eq. (A6) we can obtain the bound states. For cr = —1 
we have: 


n = 0 


= Bo{X) cosh''''‘^(a;i/(5)2Ai i, A -I- 5, 5; - sinA(a;i/5) 


= 0 


(AlO) 
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n = 1, 3, 5,... 


cosh^'''^(xi/^) sinh(a;i/(5)2i^i a„ + !,&„ + 1, |; — sinh^(a;i/5) 

= BniX) cosh^(xi/(5)2Fi a„ + i, &„ + i, i; - sinh^(xi/A) 


(All) 


n = 2,4, 6, 


J 4”^ =R„(A)cosh^+d^i/<y)2Bi 

On + ^,b 

n + i;-sinh^(a:i/(5) 


1 = zC„(A) cosh'^(xi/d) sinh(a::i/(5)2Bi 

an -1- 1,6„ -f 1, |;-sinh^(a:i/d) 


(A12) 


with: 


o-n — 2 ’ 2 


and with energies po = —avp 2 for n = 0 and: 


Po = Ml -n)+ v'^pI 


(A13) 


(AM) 


for n ^ 0 (again, to obtain the solutions for cr = 1 one just 
has to interchange the chiralities). There is a maximum 
n before is positive again: N < X. Then we see that 
the number of bound states is given by the largest integer 
less than A + 1. Regarding the constants Bn and Cn, we 
impose the solutions to be normalized and obtain for the 
chiral state n = 0: 


Bo(A) = 


(-1/4)^ 


(5B_i(A,l-2A) 


(A15) 


where Bz{a,b) is the incomplete Beta function and can 
be written as: 


with: 


='lp^^'^\xo,X2)UR.L (A18) 


Again for cr = 1 one should interchange the chiralities. 
The summation is over the bound states and the integral 
is over the continuum of extended states. 


Appendix B: Non-adiabatic correction to the CS 
term 


Let us write the CS term as: 


- - - 
^^0,0 ~ 


rA(xi) 


e^^P’'dn 


To find the non-adiabatic correction F{xi) we should im¬ 
pose the anomaly cancellation. The gauge variation of 
the edge theory (eq. @) has to be canceled by the 
gauge variation of the CS term. With this condition we 
get the differential eq. for F: 


diF{xi) = 2pl{xi) 


Bz{a,b) = f F ^(1 —t)^ ^dt (A16) There is another condition A should fulfil: 

Jo 


Finally we can decompose the fermionic field as follows: 
(for cr = -1) 


lim F{xi) = ±1 

£Ci —>-±oo 


(B2) 


(B3) 


N 


N 


^ = X! ^o,a: 2 )$^Ta:i)-f^ 4'^"Aa;o, X 2 )$i (a:i)-f 


n=0 


n—1 


+ J dk (^4'^^(a:o,a;2)$fl^(a:i)-f 4'^''^(a:o,a:2)$i ^(xi)) 

(A17) 


so that in the asymptotic limit eq. (18) is recovered. 
Then F can be computed to be: 

Fxixi) = -SBq(X) sign{xi)Q B_j^2 - A, ^j (B4) 

as long as X ^ 1/2,3/2, 5/2..., with y = cosh (a:i/(5 ) and 
S the incomplete Beta function defined in eq. (A16). For 
half integer A we get: 


Bi/2(a;i) = 4(5i3o(i) arctan (tanh(|j)) 

F'3/2 (xi) = SBqM [2arctan (tanh(|j)) -|-sec/i(^) tanh(^)] 

F 5 / 2 {xi) = ^B^d) [barctan (tanh(|i)) + sech{^){3 + 2sech{^M tanh(^)] 


(B5) 
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Appendix C: Even part of the matter contribution 
in the bulk theory 

The matter contribution to the polarization function 
in the bulk can be written as a sum of an even plus an 
odd part: 

nra« = nra«.e + nZ«.o (Cl) 

The even part is a symmetric second rank tensor. As¬ 
suming rotational invariance, the most general tensor of 
this type can be constructed as a linear combination of 
Viii'y and + p^u^, where = (l^^) 

defines the rest frame of the system (see for example^. 
Imposing transversality ^(jp) = 0) we obtain the 

general form for the even contribution: 


not appear. In that case, the even part of the polarization 
function must be proportional to — p^p'^ jp^ which 
implies Gi = —G 2 - 

To simplify the computation we will treat the static 
limit (p° = 0) and long wave limit (p = 0) separately. In 
the static limit we have: 

TtOO _ / 0 _ n'l . TjOi _ TTlO _ 0 


nL«.e = G2 (p° = 0) (<5*^ - (C4) 


nz.«,e(p) = + (Gi + G2)Pr (C2) 

with Pl° = = PjO and: 


Pi = 

|p|2 


(C3) 


and where Gi and G2 are scalar functions of and |p|. 
In the vacuum there is no preferred rest frame, so can 


The calculation of Gi and G2 has been done in the con¬ 
text of massive graphen^^. The only difference with 
the present case is that in graphene there is a mul¬ 
tiplicative factor of 4 that counts the two valleys and 
the spin degeneracy. We should note that computations 
for graphene are done with cut-off regularization, which 
actually breaks gauge invariance. However the matter 
part is finite and independent of the regularization, so 
result^ can directly be imported as long as the degen¬ 
eracy is set to 1. In the static limit we have: 


Gi(p° = 0) = 


e^0{\p\ - \m\) 




— 4m^ 


farccos - arccos (-^Jl^=)9{\p\ - 2pp) 

\ W / 7 , 21,.^|2 _|_ 4 , y,2 / ^ -^u^IpP -|- 4m^ ' 


SttuIpI \ ^ \/'y^|pP + 4m2 


(C5) 


G2(p° = 0) = 


e^e{\p\ - \m\) 


?-!rtA-44«(ipi-2pri+ 

47r 47r y v^\p\^ 




— 4m^ 


8 I I (arccos A — J - arccos A / f)^(|p|-2pF) 

StT^IpI _|_ 477 .^ 2 / V-^t;2|p|2 _|_ 477 j 2 / 


(C6) 


where the Fermi momentum is defined as vpp = 
\/— rrP- In configuration space, restoring the Xi de¬ 
pendence of m we have: vpp = p? — tanh(a;i/(5)^. 

When \p\ > Too we are safe to do 0(|p| — 2pp) = 0 for the 
description of the low energy theory as long as |yr| doesn’t 
get too close to toq- On the other hand, when \p\ < toq 
we will always have a region in space where Pp is small. 
However, this region is localized near the DW where the 


adiabaticity is lost. This means that within the adia¬ 
batic approximation we are safe to do 0(|p| — 2pp) = 0 
for any value of p except when \p\ ~ toq, in which case 
corrections proportional to 0(|p| — 2p^) would be needed. 


Then, supposing |/r| is below or sufficiently above TOq, 
we can forget about the terms proportional to 0(|p| — 2^^-) 
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and do a derivative expansion to get: 
Gi(p° = 0) = 0i\n\ - \m\) 


127r|r 


+ o{ 



(C7) 


G2(p° 


0) = d{\n\ - \m\) 



(C8) 


This terms are finite in the limit p —0. 

At the same time, in the long wavelength limit we have: 

no„Ve = 0; n^,„,, = G2(p = o)d*^ 


nLtM = n^a*M = G'i(p = o)^ (C9) 

Note that we maintain a linear dependence in p* as we 
impose spatial homogeneity in the electric and magnetic 
fields. Now we search for terms linear and quadratic in 
in Gi _2 which are finite in the limit p —^ 0, so that we 
go up to second order in a derivative expansion. With 
this procedure we get rid of non-local terms (negative 
powers in p). FrorrPSl -y^re obtain, in the low energy regime 

Pq << TO^: 


Gi(p = 0) = e{\p\ - \m\) 


I 127rw^|r 


(CIO) 


G2(p = 0) = -6»(|p|-|to|) 




m 771P 


y 127r7;^|r 

(Cll) 

where n.l.t. are non-local terms. 

Adding the vacuum and matter contributions (in both 
static and long wavelength limits) just obtained and ne¬ 
glecting the non-local terms in eqs. dcioiciiD we arrive 
to the following expression for the even part of the po¬ 
larization function (in configuration space): 


+ ,(l^l-M) d(IHr.M) . (C 12 ) 


n°» = n:° = 0(H-|p|)-^^ (C13) 

127ru|m| 


Let us note that if we do po —>■ 0 in the expressions for the 
long wavelength limit, the polarization function obtained 
does not coincide with the one we get if we do p —>■ 0 in 
the expressions for the static limit. This is due to the 
non-locality of the matter part. Withing our approxima¬ 
tions, the commutator of the limits po —>■ 0 and p —>■ 0 
gives precisely the constant term of Gi(p° = 0), which 
contributes to If™ in the static limit. This means that 
when adding the static plus the long wavelength contri¬ 
butions, non-local terms which are constant in the limit 
Po —> 0 and zero in the limit p —> 0 are being approxi¬ 
mated by a constant term. 


Appendix D: Action for the in-plane magnetization 


To obtain the action for the in-plane magnetization we 
can directly import results from section]^ The only place 
where a bit of care is needed is in the computation of the 
action for the chiral edge theory, as the sign of the chiral 
anomaly depends on the sign of the mass. For clarity we 
will obtain the first order and second order terms (in the 
magnetization) separately. 


1. First order terms 


From eqs. ( |10|12|24|25| ) we obtain: 


p(i) 


J dtdxdy Aj;y s • nixy 


(Dl) 


n 


00 

e 


-0{\m\ - IpI) 


e^\d\^ 

127r|TO| 


where the spin density s is (restoring h)\ 


s = 


pjix) 

h 


\vp OoAavpOy 


x-f 


2hvp 


aF\{x) 9(^Al tanh^(a:/d) — p" 


aAz tanh(a:/d) 

\P\ 


— Al tanh^(a:/(5)^ 


E 


(D2) 
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It is related to the electromagnetic current density as: 
j = X z. To evaluate the nonequilibrium term, 

which is non-local (the one with inverse derivatives), we 
will assume the electromagnetic current to be time inde¬ 
pendent. We gelP^ 


wall. The nonequilibrium edge current density along the 
DW then reads: 

(D5) 


± 


1 


do ± avpd, 


■Ev = - — Viy) 


VF 


(D3) 


with V(y) = —Ey y + constant. We will further set the 
voltage to zero at y = —T/2, where L is the length and 
y = EL 12 are the end points of the DW. We finally have: 


± 


do ± avpdy 


Ey = - — AV{y) 


Vf 


(D4) 


where AV{y) = V{y) — 14(—L/2) is the voltage between 
the end point y = —L/2 and a given point y along the 


For more clarity let us take the spatial component and 
compute the average 

O'L) = T / ^ = ±0-|vlPA(a:)pAF (D6) 

where here AV = V{L/2) — V{—L/2) is the voltage be¬ 
tween both end points of the DW. A similar curre nt co n- 
figuration has been used recently in the literatur^SMSI, 
Finally we can write the spin density as: 


s = (/i - eAV (y)) 


^ — [ (jFx{x) 0(AI tanh^(a;/i5) — -I- tanh^(x/(5)) | E 

2hvF \ ' \y\ ^ 




2. Second order terms 


the static limit) and (251 give corrections to the exchange 


The quadratic terms in the magnetization can be ob¬ 
tained from eqs. ( |12|24|25| . Expressions ([M| (taken in i®’ 


energy and the Berry phase term respectively. Their con¬ 
tribution to the second order terms of the effective action 


t(2) 

hulk 


= -d 


dtdxdy < tanh^(a;/i5)) (^{d^rux)^ + [dyruy)^^ 


T 


T 


Mfff 


(^aFx{xi) d(Al tanh^(a;/(J) — y? 


crA^ tanh(a;/5) 


— A1 tanh^(a;/5))^ (jnydtnix — mxdtniy'j > (D8) 


with: 


= 


A2 

xy 

127rA2C?’ 


* 2hdvl 


(D9) 


The values of the effective exchange constant and satu¬ 
ration magnetization are « 4.24 x and 

~ 4.3A/m, which is much smaller than the values 
of the ferromagnet A and Mg, so that the contribution 


h2) 

bulk 


can be neglected. Besides the exchange energy and 


saturation magnetization renormalization, eq. (24) gives 


also an extra dynamical contribution, but is second order 
in time derivatives, and hence higher order in the deriva¬ 
tive expansion than the CS contribution and can be ne¬ 
glected. There would also be further non local dynamical 
corrections at finite density, which are highly non trivial 
to compute and which should be again of no importance 
for the physics compared to the Berry phase term of the 
ferromagnet. 
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Regarding the second order contribution coming from 
the edge theory (eq. (HI) we can write: 


.(2) 

anomaly 


A? 


2hvF 


dtdt'dxdx'dydy'px{x)px{x') x 


= 2 sign{t - t')5{±avF{t -t')+y- y') (Dll) 


X nixit, X, y){-idt)G{t -t',y- y')m^{t', x', y') (DIO) 
where we defined the green function G as: 


G{t-t\y-y') = J 


47r2 go T crVFQy 


Doing the derivative of the Green function we get: 


-^(2) 

anomaly 


A2 

xy 

2hvF 


dtdxdx' dy p\{x)p\{x') m^it, x,y)mx{t, x\y)+ 


J dt'dy'm^{t, x, y)sign{t - t')dt5{±avF{t -t') + y- y')m^{t', x', y') 


r 


(D12) 


Taking into account that the magnetization does not 
depen d on t he y coordinate, the second term of the r.h.s. 
of eq. (D12) can be simplified to: 


A^ f 

±a / dtdxdx'dy p'l{x)p\{x') mx{t,x)x 

ihVF J 


X (jrixitT L/vf,x') +mx{t± L/vf,x') - 2mx{t,x')j 

(D13) 

where L is the DW length, and L/vf is the typical time 
that takes an electron to travel the whole length of the 
wall. Making use of the translation operator we can 
write: 

I L Q 

mx{t ± L/vf,x) = e *mx{t,x) = 


L 
Vf 

so that we have: 

A2 l2 

‘ —‘'Tl/ J-'' 


= (l± —dt + + 0{^df))mxit,x) (D14) 

^ TJ P AV T? U IT' ' 


L3 


±(J 




AhvF Vp 


dtdxdx'dy\p\{x)\ |/9A(2;0rx 




X mx{t,x)dtmxit,x') + 0(-Fdf) 

vp 


(D15) 


This is second and higher order in derivatives so can be 
neglected, and we finally get: 


■^(2) 

anomaly 


= J dtdxdx'dy p\{x)p\{x')x 


X mx{t,x)mx{t,x') (D16) 

where the effective hard axis anisotropy constant is: 

A1 


= -j-P- « 3.49 X lOV/m^ 
dohvF 


(D17) 


Assuming a ferromagnetic thin film with a weak in¬ 
plane anisotropy, the whole contribution to the hard axis 
anysotropy can be assumed to come from . 


Appendix E: Action for the out-of-plane 
magnetization 

Let us now compute the contributions linear and 
quadratic in m^. We integrate out the fermions in eq. 
with magnetization given by eq. (32), which in¬ 


cludes fluctuations around the equilibrium configuration 
(eq. (HI)). We work in the adiabatic approximation tak¬ 
ing the mass m = ±crAz tanh(x/(5) as a constant and 
at the end restoring the x dependence. Picking only the 
terms linear and quadratic in fhz we have, in imaginary 
time: 


^ (( 90 +Po ~ T ^^(qT(a- f i?-f C) (El) 
where £(q) = i/lqp -|- and: 

A = —AlTr — fhz{p)rhz{—p) (E2) 
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B = i/s., Tr 


(i{f + f)-m)<f{p){i^-m) rh^{-p) (E3) 


+ 2mqaap{p)5°'^ 


+ mp • a 



(E6) 


C = iA^ Tr 


{i{f + f)-m){if-m)(^{-p)m^{p) (E4) 


We defined q), p^ = and 

^ s’'5^1, since we are working in euclidean spacetime. 

It can be seen that at p = Q the calculation gives no dy¬ 
namical terms first order in time derivatives, so in this 
limit dynamical contributions can be neglected. At finite 
p we will take the static limit, assuming any dynamical 
contribution arising from finite density effects can again 
be neglected compared to the Berry phase term of the 
ferromagnet. Hence doing Po = 0 and performing the 
traces we get: 

A = A^(^2qQ+2lql^-2m^ + qtPjS"^^m:;(p)m,(-p) (E5) 


B = 2A,[e°'^^ qopaapip) + e°“^(7aP/3ao(p)-b 


C = -2A,(^e°°'^ qoPaap{-p) + qaPi5ao{-p)- 


- 2mqaai3{-p)d°‘^ - mp- &{-p)^m^{p) (E7) 


The first term of B and C has the form of a DM inter¬ 
action nixy ■ prhz, while the second term would give a 
coupling to the electric field. 

We go up only to linear order in the external momen¬ 
tum, so we expand the denominator of eq. (El) to first 


order in p and pick only up to linear terms in the re¬ 
sulting expression. Then performing the integrals in the 
internal momentum and using dimensional regularization 
for the divergent integrals present in the vacuum contri¬ 
bution we get (still in imaginary time): 


r. = 


(Pp m{p^ — m?) 


IGttuf J {2tt)^ 




^(/i^ - m?) p • {a{p)mz{-p) + a(-p)m^(p)^ -f 


A2 

Z 

Sttv], 


Pp 

(27r)3 


^A\m\e{m‘^ - p^) + (^\m\ + \p\ + - m'^)^mz(p)mzi-p) 


(E8) 


r 


The first integral of eq. (E8) above is antisymmetric in 


p and vanishes up to a total derivative. The second one 
gives the final non zero result. We see that the terms 
mixing with vanish, which means that there is 
no coupling of the out-of-plane magnetization with the 
electric field to this order and that there is not effective 
DM interaction. The reason these terms vanish resides 
in the following vanishing integral: 

/ OO 

|^((9o - ip)'^ +£'^{q))~^{qo - ip) = 0 (E9) 


Going back to real time and configuration space, and 
doing the substitution m = icr tanh(a;/i5) we get 
(restoring h): 


^/ve// ft 

=--— / dtdxdy I 4 I tanh(x/(5)| 0(A^ tanh^(a;/(5) — 


-b + I tanh(x/5)| + 


2 Az tanh^(x/^) 

P\ 


0(p^ — A^ tanh^(x/<5)) 


(ElO) 


which gives just a renormalization of the the easy axis anysotropy energy with: 

ttA? 


« 1-05 X lOV/r 




(Ell) 
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This is much smaller than the easy axis anysotropy con- stant of the ferromagnet, so this term can be neglected. 
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